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SELECT  PRINCIPAL  PARAMETERS 


DETERMINE  COORDINATES  0?  7XNGITUJ)INA 
EFFECTIVE  MATERIAL 
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INITIAL  SCANTLING 
ASSUMPTION 


BEFD3E  PIATING  THICKNESS 
BASED  UPON  LATERAL  LOAD 


INCREASE  50  m  ja  3  BR 
THICKNESS  IF  DEFICIENT 
IN  PRIMARY  STRESS 
(CYCLE  "A") 
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TSCJREASE  THICKNESSES  OF  PLATES 
SUSCEPTIBLE  TO  BUCKLING 
(CYCLE  "B"} 
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PRINT  RESULTS  OF  EACH  CYCLE  AND 
ACT  UPON  DESIGNER'S  DECISION 
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:  attempt  has  be-^n  ate  a  b 

;he  desi#aes».  general  enou&l  /et 

lie  lent  c        ^age,  Ffeasilj         the 
internal  operation  of  the  pro&raa  is  aot  aeces&r.         scanning 
the  manual,  the  user  will  gala  confidence  in  the       .deration  „ 

!E»e  seems  o:?  plating  are  numbered  in  ascending  orct 
the  keel  plate  center  up  the  shell  side;  aad  continuing  1  e  1st 

deck  to  the  inner  bottom,  if  that  exists;  and  finally  tar     he 
bulkhead  seams  o  Tlse  ship  has  been  b:     down  into  sec       nd  the 
longitudinals  within  those  sections  are  indexed  in       owing 
manner. 

TUe  bottom  shell  is  .labeled  as  sect 
longitudinal  being  indexed  as  (10.1)  and  ec  total 

number  of  longitudinals  for  toat  section.  -She  n:         up  the 
shell  side  Is  labeled  section  #11,  the  next  is  #12.     luing  until 
the  male  deck  is  reached.  !Eae  decks  are  labeled  as       fl  for 
the  main  deck,  section  #2  for  the  second,  continuing  fox     he  decks, 
the  inner  bottom  is  indexed  as  one  more  than  the  number  c     »<> 
Again  the  longitudinals  follow  the  sams  logic  as  the  rids,  l^erefore, 
(3>7)  signifies  the  seventh  longitudinal  on  the  third  deck.  Ike 
bulkheads  begin  with  the  index  of  "21"  for  the  first  bulkhead,,  aod 
continue  for  all  bulkheads.  Again  the  longitudinals  are  labeled  as 
before.  Before  operating  the  program ..  a  sketch  should  be  prepared 
so  any  confusion  will  be  alleviated. 
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1° 

The  actual  opsxation  begins  witn  the  loading  of  the 
program  into  the  computer  core  as  shows  below. 

loadgc  im  stf  piff  m?  sub  cyc  buk  mm 

STEP  2o 

Presently  the  program  Is  set  to  design  a  merchant  ship,  and 
all  principal  dimensions  are  arranged  to  fit  A.B.S.  ships.>  The  input 
data  existing  now  can  be  modified.  T3iea:?.sr0uter  will  print  a  consae 
shoving  what  symbols  can  be  changed  at  each  interval;  and  the  changes 
may  be  made. 

This  is  done  by  setting  the  symbol  representing  the  quantity 
to  be  changed  equal  to  the  new  value.  For  eyjsaaple,  (LBP  =  500)  means 
that  the  length  between  perpendiculars  will  be  500  ft.. 

Hie  information  telling  which  sections  are  to  be  transversely 
franca  is  contained  in  a  vector  called  ZRANSV.  If  sections  1,  2,   and 
11  are  to  be  transversely  teemed,   the  array  is  filled  by  typing 
,!TBAHSV(l)  ■  1,  2,   11" .  Any  section  number  may  be  included  on  the  right- 
side  of  the  equal  sign. 

Longitudinal  bulkhead  Information  is  contained  in  a  two- 
dimensional  array  called  "3HD  (I,  J)M.  The  rows  (l)  are  the  sequence 
of  bulkhead  (i.e.,  first;  second)  while  the  columns  (j),  contain 
data  describing  the  location  of  the  bulkhead.  If  the  operator  wanted 
two  bulkheads,  and  the  second  were  10  ft,  from  the  center  line,  lying 
between  decks  #1  and  #3,  he  would  type  "BBB(2,l)  =  10,  1,  3." 


An  inner  bottom  i3  provided  by  typing;  IB  =  $YES$j  if  go 
inner  bottom  is  desired  then  type;  IB  -  $H0$.  As  each  interval  of 
data  inpjt  is  finished,  the  operator  must  type  an  asterisk  (*). 
Reference  to  the  sample  session  at  the  console  will  clear  up  any 
confusion  found  in  the  manner  in  which  data  is  introduced. 

At  two  intervals  in  the  design  of  the  coordinate  system,  the 
operator  may  stop  the  program  to  look  at  the  scantlings.  2fcis  provision 
is  only  to  be  used  if  the  coordinate  system  is  being  revised,  and 
provides  a  quick  look  at  the  locations .  In  any  event,  they  will  be 
printed  out  at  the  end  of  the  program^  so  it  is  best  to  skip  them  now. 

Control  is  removed  from  the  operator  until  the  end  of  each 
cycle,  whereupon  he  can  make  his  decision  as  to  which  set  of  thicknesses 
to  use  for  the  next  cycle;  or  if  the  design  is  satisfactory,  be  may 
choose  the  final  thicknesses  and  stop*  The   final  scantling  design  is 
printed  end  the  program  is  ended. 
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SM£LE  SESSION  at  console 


loridno  mn   pit   stf    tap  eye   huk   ond    sub 
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EXECUTION, 


THE  VALUES  BELOW  MAY  HE  ALTERKO  , 


LBP,R,n,H,NnKS,HEEL,nPRiSE,FRSP,WERSP,SI      F,NBH 
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lbp*550/h=7.S/d«39.3<,h=27.5,cx«=.98/eb«.fi7^frsp«30,w?^so=in,ndks: 
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transv(l)«l,2,3,«*,10,ll,12* 


THE  VALUES  BELOW  MAY  RE  ALTEREO  NOW 
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A  decision  branch 
( "whenever  statement ) 


J <A,b  and 
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An  iteration;  for  example,  the  "thro'jigh  loop' 
where  I  is  initially  set  to  1  and  the  itera- 
tion is  performed  until  I  is  greater  than 
TOTSTRo 


The  enclosed  values  are  calculated <. 


KJK. 
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Subroutine  (BUK)  is  executed. 
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Ccatimmfcicm  from  previous  page 


Continuation  go  next  page 


Print  out  of  data 


L 


Read  :Ln  Data 


CZD 


Beginning  or  end  of  each  subroutine 
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(MN.)     s  input  and  calling  program  that  executes  the 
subroutines  and  contains  all  of  the  necessary  input  data  for  the 
operation  of  the  program., 

(MS.)  first  clears  ail  simulated  tapes  of  any  existing  data 
computed  from  previous  runs  and  theu  sets  the  values  for  the  following: 

Young's  modulus!  length  'between  perpendiculars ,  bean,  depth, 
draft,  number  of  decks,  heel  angle,  deadrise,  frtuae  spacing,  web 
spacing,  midship  section  coefficient,  height  of  center  vertical  keel, 
block  coefficient,  yield  stress,  number  of  bulkheads,  the  depth,  areav 
heigh ,  of  neutral  axis,  moment  of  inertia  and  section  modulus  of  the 
ba.3lc  stiffener,  bending  rooLoeut,  limiting  bending  stress,  bilge  radius, 
number  of  longitudinal  bulkheads  in  the  midship  half  section,  and  the 
width  of  cargo  hatch. 

The  required  ship  section  modulus  is  determined  according  to 
the  Load  Line  Regulations  and  is  calculated  by  the  subroutine  (REQSEC.) 

To  change  any  of  the:  above  existing  daoa,  there  are  two  program 
stops  which  permit  new  input  data  to  be  read  into  the  computer.  Refer 
to  the  section  on  program  operation  for  further  clarification  of  the 
main  program. 
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urn 


ABASIC 
B 

3ILL 

BM 

CB 

C3ILGBS 

CBR 

CVK 

cx 

D 

DBASIC 
DDRISE 
E 

EM 

FL30T 

FRSP 

H 

HATCHO 

HB 

HEEL 

IB 

IBA3IC 


Array  of  numbers  from  (2),  Sect,  ^3.15  -  17(c)(5) 
used  in  calculating  ZEE- 

Aspea  of  basic  stiffened. 

Breadth  of  ship. 

Linear  length  of  bilge. 

Bending  moment  of  ship  on  wave. 

Block  coefficient. 

Height  of  bilge  center  above  keelo 

Standard  camber. 

Height  of  inner  bottom. 

Midship  coefficient. 

Depth  of  ship. 

Depth  of  basic  atiffener. 

Dead  rise  of  ship  hull. 

Young's  modulus. 

Function  used  in  a  buckling  formula  later  in  the 
program. 

Length  of  flat  portion  of  bottom  shell. 

Nominal  longitudinal  spacing. 

Draft  of  ship- 

Half  breadth  of  hatch  opening- 

Half  breadth  of  ship. 

Heel  of  ship  used  in  calculating  inclined  water  head= 

Inner  bottom 

Inertia  of  basic  stlf fener 
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LIST  0?  SYMBOLS  C  : 


LBP  Length  between  perpendiculars. 

;$BEDS  Number  of  longitudinal  bulkheads  on   half  section. 

NBX3  Number  of  strength  decks  o&  ship. 

PHI  Angle  of  deadrise. 

3.I&L6 

H  Bilge  radius. 

SIGrLIM  Limiting  primary  bending  stress. 

SIQ1AY  Yield  strength  of  steel  used. 

2RANSV(I)    Array  consisting  of  the  numbers  of  sections  to  be 
transversely  framed, 

TWD         'Tween  deck  height. 

MSBSF       Spacing  of  web  frames. 

YBA3IC      Distance  from  shell  to  neutral  axis  of  basic 
stiffener. 

ZREQ        Section  modulus  of  ship  required  by  (l)o 
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EaP. 


(TAP.)  is  a  subroutine  wilier  reads  and  writes  FROPCO  (I, J), 
and  XLL(I,J)  onto  simulated  tapes  at  various  points  in  the  program  in 
order  to  conserve  storage  locations  In  the  memory  car?)  of  the  computer. 


BEQSEC. 

(REQSBC.)  is  an  internal  function  which  calculates  the 
required  midship  section  modulus  according  to  the  Load  Line  Regulations. 
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WOTET  TAPE 
VC    FR<J>M 
ARRAY  X, 

(R-fcAliNG, 


write:  tape 


array  x, 
4>fc  Y,  (J>K  %. 


fcerrue.M 


READ  TAPS 
ARRAY  Kj 
(FUWSTlMSj 


*tV  f^4>m 

ARRAY    X, 
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HEAD, 


(BEAD.)  is  an  internal  function  which  calculates  water  head, 
using  the  horizontal  and  vertical  coordinates  of  a  point  w   the  hull 
geometry  as  inputs. 

Th«  internal  function  ia  executed  by  (STF.)  and  {PLT.)>  «*"* 
calculates  two  values  cf  head  for  each  coordinate:  fh«>  head  when  the 
ship  is  heeled  thirty  degrees  and  the  head  in  way  of  the  wave  crest. 
The  larger  of  the  two  values  is  selected  as  the  water  head  and  it  in 
not  to  he  less  than  5 . £5  feet . 

L3ST  OF  SYMBOLS  CALCUIATED  IS  (HSAS. ) 

IKC1EP    Water  head  on  a  structural  mmber  when  the  ship  is  heeled 
30  degrees 

WHD      Wo/;er  head  on  a  structural  member  in  way  of  *.  wave  crest, 

HM1B      Minimum  water  head  of  5.2*  feet. 
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STP. 


{SEP. )  calculates  the  longitudinal  trm&  spacing  in  each 
section  of  the  ship,  determines  the  horizontal  and  vertical  coordinates, 
measured  from  the  keel,  of  the  longitudinal  stiffeaers,  aud  calculates 
the  water  head  an  each  stiffened  locatisa. 

The  girths  of  each  section  are  first  computed  and  then  a 
noaisal  longitudinal  frame  spacing  3s  deterained  so  that  the  longitudinal 
stiffened  are  equality  spaced  an  each  section  and  as  near  as  possible, 
hut  not  exceeding,  the  input  frame  spacing. 

Hext,  the  longitudinals  are  located  on  each  section  based  en 
the  spacing  for  that  section,  and  the  XLL  {horizontal)  and  ILL  {vertical  j 
coordinates  are  calculated.  With  these  coordinates,  the  water  heads  on 
each  stiff ener  are  computed  by  the  subroutine  (HSA2).). 

Por  sections  that  are  transversely  framed,  thus  having  no 
longitudinals,  Xli,  YLL  and  the  water  heads  for  toe  sections  are  set 
equal  to  zero  at  the  end  of  this  subroutine. 
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SYMBOLS  CALCULATED  BiffeiTF,) 


ALPHA       Angle  measuring  distance  of  inner  bottom  from  center 
point  of  bilge  turn. 

BIffST       Incremented  sum  of  frame  spaces  on  the  bulkheads, 

BILffTH      Incremented  Srjm  of  frame  spaces  on  bul^e  turn. 

C^URT       3fce  accumulative  sum  of  longitud5.u&Iij  on  shell  side 

D3LTA       Incremented  angle  between  longitudinal  upper  bilge 
turn. 

DKHT(I;      Array  of  the  deck  heights  above  the  keelo 

FRDIV       The  first  estijaate  of  the  number  of  longitudinals 
in  a  section. 

GAMMk(l-J)   Angle  of  longitudinals  on  bilge  turn  up  to  inc.* 
bottom. 

GORTH(I)  Array  of  all  the  section  lengths . 

HMAX  Maxiiaum  water  head  chosen. 

HMBS  Minimum  water  head  possible. 

EBHB  Inner  bottom  half  breadch 

IMCLHD  Head  on  longitudinals  when  ship  is  inclined. 

LC<fUKT(l)  Array  counting  the  longitudinals  in  each  section. 

LHD(I,J)  The  water  head  on  each  longitudinal* 

LL10  Accumulated  sum  of  the  longitudinals  in  Section  10 . 

LL11  Accumulated  sum  of  the  longitudinals  in  Section  11. 

LLDK  Accumulated  sus  of  the  longitudinal*;  on  the  decks. 

LL1B  Accumulated  sun  of  the  longitudinals  on  the  inner  bottom  > 
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SYMBOLS  (2A1CULATED  IN/STff. 


2  • 


S(l)        Array  of  all  the  frame  spacings  in  each  section. 

SDIST       The  incremented  frame  spacing  of  each  section  on 
the  shell  side, 

SDVAR       An  integer  which  has  the  value  NDKS  »&en  the  ship  has 
no  inn«r  bottom  and  the  value  NDK3+1  when  an  inner 
bottom  exists. 

SIDE        A*:  incremented  length  at  bilge  bottom,  stellar  to 
BILffTHc 

SUXffH      An  incremented  length  in  section  11  on  the  bilge  turn, 

SUM        A  quanity  describing  several  incremented  totals  used 
through  out  the  program. 

TDIST       Incremented  frame  spacing  seen  on  the  top  decko 

TS1PL       The  distance  between  the  last  longitudinal  on  the 
bottom  shell  and  FLBQT* 

'StyOALL  The  total  number  of  longitudinals  on  the  half  section , 

WHD  The  wave  head  on  the  longitudinals  * 

XLL(1£T)  "X"  coordinate  of  longitudinals. 

YLL(I,J)  "Y"  coordinate  of  longitudinals. 
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"PUS"   is  a  subroutine  which  determines      3  geometric 
properties  that  depend  upon  the  plate  seaias»  When  execution  of  this 
subroutine  has  been  completed;  the  coordinate  system  for  ail  losgi- 
tuditial  1 ratinuous  material  i&  the  saidciiip  section  is  fixed  and  the 
program  is  ready  to  begin  the  scantling  selection  procedure. 

The  first  step  iz\  the  execution  of  the  "PLT '  is  the  fixing 
of  plate  widths.  The  keel,  center.,  margin.,  sheer,  top.  and  stringer 
stroke  widths  have  been  calculated  in  Ref  >  (  )  according  to  ABS  rules. 
The  relations  obtained  are  v.sed  by  this  subroutine  iio  2>r3set  the  values 
of  these  plates.  Next,  a  standard,  plate  width  base^     ship  length 
is  selected  for  all  remaining  pj^ting. 

When  a  standard  size  has  been  selected,  the  next;  step  is 
test  all  seam  locations  for  interference  with  the  stif.teners  whose 
coordinates  have  been  fixed,  The  keel  plate  la  checked  first  and 
corrected  if  interference  occurs,  then  the  center  stroke  is  made  to 
equal  the  keel  width.  The  remaining  standard  plate  widths  must  now  be 
changed  to  a  nominal  size  which  yields  an  integer  number  of  plates  on 
the  shell,  the  decks  and  on  the  bulkheads.  Thus,  e     .trske  width 
has  been  determined.  It  remains  now  to  cheek  each  one  for  inter- 
ference with  stiffeners  and  to  store  aT^y  all  data  needed  by  subse- 
quent subroutines, 

The  shell  seams  are  first  examined.  The  distance  along  the 
shell  from  the  center  of  the  keel  for  each  strake  top  is  found;  checked 
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against  stiff ener  interference,  and  changed  if  necessary.  Next  the 
coordinates  of  the  seams  and  the  water  head  on  thein  is  calculated. 
The  area,  moment  and  inertia  constants  are  also  calculated  for  use 
later  in  the  program . 

Deck.,  inner  bottom,  and  bulkhead  plates  are  '.ixaoiaed  in  a 
sequence  similar  to  that  of  the  shell.  The  distance  of  each  seam  from 
the  deck  side  or  bhd.  bottom  is  found  and  checked  against  interference. 
Then  the  coordinates,  water  head,  and  area  moment  and  inertia  constants 
are  found*  Once  these  quantities  are  calculated,  the  program  returns 
control  to  the  main  program.  Before  this  is  done,  hoxv*ever,  it  is 
possible  to  obtain  a  print  out  of  the  data  calculated  in  this 
subroutine.  Refer  to  the  chapter  on  program  ejcecution  for  the  correct 
method. 
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BDSSR(I) 

C0NT1 

C'IRSTR 

ekstrCi) 

E2A(j) 

ESEAM(J) 

rc(J) 

keel 

Margin 

MC(J) 

REMPLT 

SEAM(J) 

SHEER 

SHELL 

SHL3TR 

3PANST 

STDPLT 

STRSTR 
SUM1 
TODK 
TOTSTR 


Area  eons        J"th  plates  depended  on       ear 
Length  of  j  late* 

Number  of  ctrakes  on  "I"th  bhd. 

Key  to  print  additional  ds 

Width  of  center  strake  on  inner  bottom 

Number  of  strakes  on  sI"th  deck. 

Angle  between  seara  aa  bilge,  similar  to  "GAMMA"  in  "STF" 

Water  bead  an.   "J"ta  seam 

Inertia  coaa3taat  of  "J"th  plate 

£eel  plate  width 

Width  of  margin  plate  on  inner  bottom 

Moment  constant  of  :,J"th  plate 

Length  of  ell  plating  on  shell  except  ICES.,  SHEER,  TPSTR. 

Itoe  linear  distance  from  the  top  seam  oi       plate 
to  the  Keelo 

Width  of  sheer  plate. 

Linear  length  of  all  plating  on  shell 

Number  of  plates  on  shell. 

Bbe  plate  to  be  associated  with  each  stiffened. 

A   nominal  assumed  standard  width  of  plate  to  be  dependent 
on  ship  length. 

Stringer  plate 

Accumulated  sum  of  GURTH(l)  up  the  side  from  the  keel. 

Number  of  plates  on  all  the  decks 

Accumulated  sum  of  all  the  plates  on  the  half  section > 
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itof  S^l:c  3  Calcy. 

TPSTR  Stop  strake 

WBDSTR(I)  Bominal  width  of  each  plate  on  xI"th  bhd. 

WQKSTR(l)  Nominal  width  of  each  plat 

WSHLST  Nominal  width  of  each  plate  on  st 

X  SEAM  (J)  "X"  coordinate  of  ,:J"th  seam 

Y  SEAM  (J)  "I"  coordinate  of  "J"ta  seam 
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AM1STP. 

(AMISEF. )  calculates  the  total  area,  saoaent  and  inertia 
of  the  stiff ener*  about  the  base  line. 

UBiug  the  value  of  proportionality  constant  for  each 
stiff ener  loc&tioa  determined  in  {SEEKER.),  the  area,  moment  and 
inertia  of  each  stiff ener  is  calculated.  The  resulting  values  for 
each  stiffened  can  be  added,  Bine©  they  were  tahcn  ©beat  the  same 
base  line,  thud  giving;  the  total  area,  moment  and  inertia  of  all  the 
stiff enera . 

LIST  OF  SYMBOLS  CAZ£UXA'EI£D  IS  (A&LSgP. ) 


ASEF     Tot*l  cross*sectional  area  of  the  loagltudiral  stlff#aers 
at  the  midship  section. 

2STF     Moment  of  inertii  of  the  cross s -sectional  area  ©f  the 
longitudinal  stiffeney*  about  the  base  line. 

MSTF     Internal  bending  ooEent  of  the  crcss-sectiosal  area  of  the 
longitudinal  stiff  eners  about  the  base  line. 
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(AM3PLT. )  calculates  the  total  area,  moment  and  inertia  of 
the  plates  about  the  base  lice. 

Using  the  values  of  area  constant,  moment  constant  and  inertia 
constant  calculated  in  (PIX.)»  and  the  current  value  of  plating 
thicknesses  In  the  sub=cycle,  the  area,  moment  and  Inertia  of  each  plate 
about  the  base  line  is  calculated.  The  resulting  values  are  added, 
giving  the  total  area,  moment  eatf  inertia  of  sOi  plates  about  the 
base  line. 


LIS!*  OF  SYMBOLS  CAICUIAT&i)  Ifl  (AMXPLT.) 

APLT     Total  cross-sectional  area  of  the  plating  at  the  midship 
section. 

IPLT     Moment  of  Inertia  of  the  cross^sectlonal  area  ©f  the 
plating  about  the  base  line. 

MPLT     Internal  or  resisting  bending  moment  of  the  crass -sectional 
area  of  the  plating  about  the  base  line. 
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tioa  Modulus  versv 
*b>      utine  (SI 

by  Lt&raogian  lateri 
proportionality  coaata&t  of  the  basj, 
required  section  modulus     ?late-s     ier  ooi. 

LIST  OF  SYMBOLS  CAIOU^^^^IgmP^ 
TiMP       Ifcisajy  variable  equal  to  the  proportionality  constant. 
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For 
sise  o:  »enes.  c 

tion  to  withstand  the  maximum  stress-. 

Whenever  (STFNER.)  is  executed  by  (CYC,        ordiaatea 
the  stiffener 3,  water  heads,  plate  thicknesses  assoeii     ith  each 
stiffener  location,  and  the  stresses  for  the  current  s 
It  is  the  purpose  of  this  subroutine  to  determine  the 
scantlings • 

A  scale  factor  or  proportionality  constat         :ased  on 
the  geometry  of  the  basic  stiffener  and  an  effective  plate  width 
60  x  thickness  or  the  longitudinal  t  spacing,  w       is  less, 

used  to  determine  the  geometric  properties  of  the  required  stiffener. 

For  each  stiffener  location  and  its  assoc     plate 

thickness,  (STFHER.)  calculates  fifteen  values  of  section  modulus  for 

each  plate-stiffener  combination,  by  using  fifteen  different  values  of 

PROFCC*  In  effect,  a  curve  of  section  modulus  versus  proportionality 

constant  is  generated  for  each  stiffener  location.  The  actual  section 

modulus  required  to  resist  the  stresses  in  the  current  sub-cycle  is  a 

point  on  the  curve.  To  design  a  ship  that  approximately  meets  the 

American  Bureau  of  Shipping  Rules,  a  section  moduli     tired  by  the 

Rules,  MOEABS,  is  initially  selected  as  the  minimum  recuired  sect:* 

modulus  of  the  plate-stiffener  combination.  With  this  initial  value  of 

MOBABS,  the  subroutine  (ISTERP. )  interpolates  the  section  modulus 

proportionality  constant  curve  to  determine  the  PKOPCO  associated  wit 

the  initially  selected  value  of  MOPABS.  With  this  value  of  PROPCO, 

(STFBER.)  calculates  the  properties  of  the  stiff eaer-plate  combination 

for  that  stiffener  location. 
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To  determine  if  this  combination  will  fail  due  to  ii 
under  the  combined  prims       ^ondary  stress*:     vsheako's 
Fosmila  for  combined  lateral  and  axial  Ioe  used. 

The  section  modulus,  2?1M,  from  Tisuosheniio ' s  Formula  is  a  function 
the  gecacetric  properties  of  the  plate-stiffener  combination  and  the 
applied  loads.  Thus,  corresponding  to  the  initial  value  of  MOBABS  and 
PROPCO,  a  value  of  section  modulus,  ZP1M,  is  calculated  by  Tlmoshenkc 
Formula.  Three  possible  situations  exi3t  at  this  point: 

I.  2T1M  equals  M02A3S 
II  o  ZT1M  less  than  MOBABS 
III.  2S1M  greater  than  MOBABS 

In  case  I,  when  2T1M  and  MOBABS  differ  by  no;  more  than  C.75, 
MOBABS  is  considered  to  equal  ZT1M,  and  is  the  required  section  modulus 
for  the  plate-stiff ener  combination. 

In  case  II,  since  ZT1M  is  less  than  the  ABS  value,  MOBABS 
governs  and  the  value  of  PROPCO  associated  with  MOBABS  is  the  required 
one. 

In  Case  III,  when  ZT1M  i3  greater  than  MOBABS,  ZT1M  must  be 
the  governing  section  modulus.  But  the  value  of  PROPCO  was  determined 
from  the  initial  value  of  MOBABS  and  does  not  correspond  to  ZT1M.  A 
new  value  of  MOBABS  is  chosen  such  that  the  value  of  PROPCO  determined 
by  this  new  value  of  MOBABS  also  corresponds  to  the  value  of  ZT1M.  To 
converge  on  a  solution s   the  mean  value  between  ZT1M  and  MOBABS  is 
selected  as  the  new  initial  value  of  I«50BAB3o  A  greater  value  of  PROPCO 
results,  and  the  two  section  moduli,  2&'1M  and  MOBABS  based  on  this  new 
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stiffens  scantling      sati 

procedure  is  repeated  for  each  at     r  location. 
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THKNES. 

thicknesses  of  plating  required  to  withstand  the  ata        head. 
Each  plate  has  two  possible  geometric  arrangements.  The  geometry  of 
the  transversely  framed  and  longitudinally  framed  sections  each  require 
a  separate  stress  distribution.  The  subroutine  (MAXSTR.)  is  called  to 
calculate  the  maximum  stress  associated  with  each  section;  this  will 
be  the  stress  upon  the  plates  in  the  transversely  framed  sections. 
Otherwise ,  the  stresses  calculated  in  (STBES1.)  are  used  as  a  basis 
for  the  longitudinally  framed  sections.  Those  stresses  i&  (STRESlo) 
are  associated  with  each  seam.  The  stress  on  each  plate  is  deter- 
mined from  the  maximum  of  the  two  seam  stresses. 

The  array  CO.   is  next  filled  with  the  plating  thicknesses 
based  upon  formulae  that  approximate  ABS  rule  requirements  (Ref  .  Hoc  ), 
using  the  stresses  described  above.  These  thicknesses  are  those 
necessary  to  meet  the  stress  schedule  set  up  in  the  current  sub*cycle 
being  run  in  subroutine  (CYC). 

Several  different  values  are  used  to  calculate  the 
thicknesses  of  each  plate.  Since  the  program  is  currently  set  to  meet 
merchant  criteria,  some  plate  scantlings  depend  upon  the  location,  and 
upon  the  ABS  rule  requirements.  The  formulas  for  the  center  plate, 
keel  plate,  sheers trake,  stringer  plate,  and  top  strake  are  presently 
all  special  cases  designed  to  meet  ABS  requirements.  A  re-orientation 
of  the  long  side  in  the  longitudinally  framed  section  requires  a  new 
formula  for  calculating  CO^   This  is  triggered  by  the  urogram  when  a 
section  having  no  longitudinals  is  encountered. 
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bhe  new 
thickse  ses  calculated  in  the  last  sub-cycle 

T12  =  ^COi  +  OT  x  V'  C«T  +  1)  {  1  ) 

^12  =  Hev  *lilicl£ness  of  itJl  plate- 

Til  =  0M  tllickness  of  ith  Plate. 

CO.  =  Thickness  of  ith  plate  just  calculated  baaed 
upon  local  lateral  load. 

WT  =  Weighting  fraction  associated  with  1th  plate. 

The  weighting  function  is  used  to  prevent  the  change  in 
thickness  during  each  sub-cycle  from  becoming  too  rapid,  thus  upsetting 
the  cycling  procedure. 

The  iteration  for  each  plate  in  each  sub-cycle,  once  the  CO- 
array  is  completed,  consists  of  three  3teps.  First,  a  correct  weight** 
ing  function  is  determined.  This  function  is  designed  to  permit  the 
most  rapid  change  of  plate  thicknesses  in  each  sub-cycle  to  ©ccur  on 
plates  which  will  not  be  changed  by  the  linear  extrapolation  program* 
These  are  the  plates  furthest  from  the  fibers  which  are  found  to  have 
plates  whose  primary  stress  is  above  the  limiting  steps.  The  second 
step  is  to  fill  the  second  column  (T.-)  of  that  six  column  array 
described  in  (CYC).  This  is  done  in  accordance  with  formula  (  j  )q 
At  the  first  calling  of  this  subroutine  by  (CYC, )  the  initial  thicknesses 
are  calculated.  This  is  signified  by  setting  the  value  of  COUMT  equal 
to  10.  For  this  case  there  is  no  set  of  thicknesses  in  column  one  of 
the  array  to  be  averaged,  so  column  two  is  set  directly  equal  to  the 
thicknesses  based  upon  local  lateral  loado  The  last  step  is  the  replace- 
ment of  column  one  by  the  new  thicknesses  In  column  two.  These  are  now 
the  new  initial  thicknesses  to  be  used  in  the  next  sub«cycle.  Control 

is  returned  to  (CYC.)o 
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LIST  OF  SYMBOLS  CALCULATED  IK  (THKNES.) 


CO(J) 

T(I,  CYCLE) 

2(1,  CYCLE  +  1) 


WT 
X 


The  plate  thickness  baaed  upon  lateral  load. 

The  thickness  used  toy  each  subcyele  of  (CYC}. 

The  thickness  found  from  the  weighted  average  of 

?  (I,  CYCLE)  and  CO(J), 

The  weighting  function  for  each  plate. 

The  maximum  seam  stress  associated  *ith  each  plate 
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-iave  't. 
program  must  now  d 

thj     es,  arbitrarilj 

stress;  the  scantlings  are  , 

which  meets  all  stress  criteria  in 

cycling  procedure  is  rather  in  understand  the  des.v 

in  its  entirety,  an  understanding  of  each 

the  overall  scheme  of  the  main  cycle  is  mandatory. 

The  ultimate  aim  of  each  main  cycle  in  (CTC.)  is  to  satisfy 
the  limiting  primary  bending  stress  on  the  hull  girder.  Cycle  #1 
begins  with  the  selection  of  an  arbitrary  se»  of  initial  thicknesses  as 
mentioned  above.  Tfei3  -job  is  done  by  subroutine  (THKNE3.)  based  on  a 
stress  schedule  set  up  by  (STKESl.).  The  stiffened  scantlings  are  now 
determined  by  (STFNER.)  based  upon  these  stresses  found  in  (STKES2.)° 
Subroutine  (AM1PI3?.)  and  (AM1STR)  determine  the  area,  moment,  and 
inertia  of  this  initial  set  of  scantlings  aad  now  Cycle  #1  moves  to  a 
sub-cycle  which  designs  the  scantlings  to  resist  1  ;  due  to 

static  water  head.  (STKESl.)  and  ( 3XRE32. )  set  up  a  stress  schedule 
based  upon  areas,  moments,  and  inertias  of  previously 
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asses 
to  resist  later :  ^TSR.}        .>s  the  as3ociat.; 

refinement  is  accomplished. 

The  next  step  is  to  attempt  to  satisfy  the  primary  bending 
stress  using  the  thicknesses  found  by  the  static  head  requirement. 
Shis  will  be  called  Cycle  A,  and  refers  to  the  part  of  the  cycle  which 
increases  the  thicknesses  of  plates  that  lie  on  the  cuter  fiber.  Cycle  B 
will  refer  to  the  cycle  which  attempts  to  satisfy  the  primary  bending 
stress  by  increasing  plates  which  would  buckle  under  the  current 
stress  schedule.  A  trigger  is  set  that  tells  (CYC.)  which  cycle  is 
being  executed.  Whenever  the  quanity  METER  is  set  equal  to  1,  Cycle  A 
is  being  executed,  but  when  METER  equals  0,  the  subroutine  is  perform- 
ing Cycle  Bo  Cycles  A  and  B  both  make  use  of  a  linear  program  to 
increase  the  area  of  the  midship  cros3  section  to  reduce  primary  bending. 

The  array  T.  .  i»  *hovn  In  Pig.  as  having  six  columns,  and  a 
number  of  rows  equal  to  the  number  of  plates  on  the  cross-section.  The 
first  column  is  the  thicknesses  of  plating  needed  to  resist  lateral  load. 
The  remaining  five  columns  are  set  equal  to  the  first,  in  order  that  the 
linear  prograimmay  operate  upon  the  array  to  chante  the  thicknesses  for 
an  optimum  stress  decrement. 
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T   *  Platea  bel©v  aeutral  axi»  iacreated  by  liaear 
extrapolation. 

T   =  Plates  above  neutral  axis  Increased  by  linear 
tX      extrapolation. 


L   =  Column  index 
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*fie  first  step  iu  the  linear  program  1.8  t>o  jiulc'-  e.<i 
incremented  increase        sicknesses  of  the  bottom  placing  in 
column  2  and  to  the  top  plating  in  i  te  lin«     ogram 

calculates  the  unit  effect  of  these  increments  and     .ys  the  results 
in  columns  4  through  6.  Column  k   is  the  set  of  new  top  thicknesses 
causing  a  stress  decrease  in  the  top.  Column  5  is  the  new  bottom 
thicknesses  causing  a  stress  decrease  in  the  bottom  fibers.  Column  6 
has  new  top  and  bottom  thicknesses  for  a  decrease  in  stress  at  both 
the  bottom  and  top. 

Column  k  through  6  are  checked  against  the  existing  stresses 
from  Column  1. 

(TOPSTR(l)  and  ;boTSTP.{1))  to  find  the  correct  column  to  use 
as  the  one  which  lowers  the  stress  for  the  best  attempt  to  satisfy  the 
limit  stress. 

If  Cycle  A  is  being  run,  the  top  (Tt,)  or  bottom  (T.,) 
sections  of  the  array  are  plates  which  are  specifically  designated. 
They  are  the  bottom  plating  up  to  the  bilge  top,  the  shear,  top  and 
stringer  strakes,  and  the  main  deck  plating.  In  Cycle  3,  the  plates 
are  the  top  plates  which  would  buckle  and  the  bottom  plates  which  would 
buckle.  During  Cycle  3,  subroutine  (BUK.)  i;  called  to  find  which 
plates,  if  any,  would  buckle  for  use  in  the  linear  program.  The  secorjL 
time  it  is  called  to  handle  a  case  when  the  limiting  s ;resa  is  satis- 
fied, but  some  plates  still  buckle. 
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Since  Cycle  A  dent,  thfc  inforj&atiou 

gained,  from  Cycle  A  aust  be  stored  wails  Cycle  3  is  run.  Upon  complex 
of  Cycle  B,  tiae  results  of  both  A  and  B  are  printed  for  the  designer  to 
make  a  decision  on  which  set  of  thicknesses  should  be  ueed  c.a   the 
initial  set  in  Cycle  2.  If  the  deal®     sss  is  to  be  terminated, 
subroutine  (ETCh)  is  called  fo?r  final      out  of  the  design.  For 
recycling-  all  necessary  quanitie3  are  initialized  and  control  is 
transferred  back  to  the  refining  subcyeie  for  further  lateral  instability 
calculation. 
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EA 
AREAI 
AK3A2 
AR(I) 
ASTIF 

snifi 

HDIf2 

BOTSTE(l) 

COUHT 

CYCLE 

CYCLES 

DEF1 

DEF2 

BET 

DT 

INEKT1 
INERT2 
IBRTIA 

WD 

ISTIF 


al  area  result!? 

Area  of  plates  in  each  c<  ,J). 

Temporary      je  of  &     ,or  arc 

Maximum  differ       een  re<  al 

plate  buciLling  stress  In  cycle 

MaxlnuE>.  differeac  ctual 

plate  buftkling  stress in  cycle  . 

Beading  stress  on  bottom  fibex       ompletion  of 
the  refining  sub-cycle . 

A  counting  device  to  trigger  certain  operations. 

A  count  of  the  desired  column  1     J). 

The  current  man  cycle  of  iters 

Current  deficiency  in  top  stress 

Current  deficiency  in  bottc 

A  determinant  in  this  linear  pi 

Small  increment       kness  of  plates  at 
designated  locations  for  linear  program. 

Total  moment  of  inertia  at  end  of  cycle  A. 

Total  moment  of  inertia  at  end  of  cycle  B. 

Current  tonal  moment  of  inertii     cross  section. 

Moment  of  inertia  of  each  column  in  T(I,J). 

Temporary  storage  of  stiff ener  inertias. 
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1MDXFQ 

METER 

MODDK 

moddki 

MODDKL 

MODKL1 

MOMENT 

MO(I) 

MOM1 

MQM2 

SlGSl(l) 

SlG13V(l) 

SlGL3V(l) 

STOP1 

TL1M(I) 

TOPSTB(l) 

TSAVE(I) 
TS(I) 

Xll 

*21 
*22 


Deficiency  in  primary  beading  stress  at  end  of  eyclf 
Deficiency  in  primary  bending  stress  at  end  of  cycle  B. 
A  counting  device  showing  whether  cycle  A  or  B 

Section  modulus  to  deck  in  cycle  A. 

Section  modulus  to  deck  in  cycle  B. 

Section  modulus  to  keel  in  cycle  A. 

Section  modulus  to  keel  in  cycle  B. 

Current  total  moment  of  cross  section. 

Total  moment  of  each  column  in  T(I,J). 

Total  moment  at  end  of  cycle  A. 

Total  moment  at  enci  of  cycle  BB 

Stress  on  each  plate  from  cycle  B. 

Stress  on  each  plate  from  cycle  A. 

Initial  stress  on  each  plate  in  each  main  cycle. 

Trigger  to  act  upon  designer's  decision. 

The  column  of  T(l,J)  used  in  cycle  A. 

The  bending  stress  on  the  top  fiber  upon  completion 
of  the  refining  sub-cycle3. 

The  initial  thicknesses  for  cycle 

The  initial  thicknesses  saved  for  cycle  Bo 


Coefficients  of  the  restraint  matrix  in  the  linear 
program. 
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coluinn  4, 

Y2{4)  Actual  inc 

in  column 

^1(5 5  Actual  increase  in  thickness  of  the  top  fibers  in 

eoluxan  5» 

y2(5)  Actual  increase  in  thickness  of  the  bottom  fibers 

*ai  column  5. 

yl(6)  Actual  increase  In  thickness     e  top  fibe 

column  6. 

Y2(6)  Actual  increase  in  thickness  of  the  bottom  fibers 

in  co?.' 

Y8A  The  neutral  axis  of  the  cross      .  \  at  ^he  end 

of  each  c,vcle. 
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Thi*  subroutine  tents  the  plate  thicknesses,  obtained  from 
local  lateral  leading  *cr  :La$tability  under  the  applied  compressive 
stresses. 

The  subroutine  Is  called  by  (CYC, )  and  consists  of  two 
mala  parts:  one  that  tests  the  longitudinally  framed  sections  far 
buckling  and  tbe  ether  tbat  tests  the  transversely  framed  sections 
for  buckling. 

For  sections  that  ara  longitudinally  framed,  the  critical 
buckling  stress,  SXGBUL,  ie  determined  for  each  panel  by  Bryan's 
Formula.  The  current  panel  thicknesses  are  used  in  the  fonsulr,  and 
when  two  plate  seams  lie  between  longitudinals,  the  subroutine 
selects  the  thickness  of  tlie  thinner  plate  to  be  used* 

The  actual  primary  stress,  31GBHL,  is  calculated  as  the 
stress  on  the  point  of  each  panel  furthest  from  the  neutral  axis. 
When  S1GBKL  is  less  than  S1GBUL.  the  panel  thickness  remains  as  the 
current  value  TB.  When  S16BKL  is  greater  th*a  S1GMJL,  the  current 
panel  thickness  is  increased  to  the  mean  value  between  the  current 
panel  thickness,  TB,  and  the  required  thickness,  TL(J),  the  thickness 
required  to  withstand  SIGBKLo  When  S1GEBKL  has  been  decreased  to 
within  100  p<,8.i„  of  SiGBUL,  the  panel  is  considered  to  satisfy  the 
buckling  stress.  The  factor  of  100  is  used  to  prevent  excessive  cycling 
in  the  computer. 
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For  sections  that  are  transversely  framed,  the  primary  stress 
acting  e&  a  panel  is  calculated  by  the  subroutine  (MAXSTR*)*  which 
determines  the  maximum  compressive  stress,  S1GEK3?,  acting  gs  each 
transversely  framed  panel.  The  critical  buckling  stress,  S1GBUT,  is 
calculated  by  Montgomerie •  s  Formula*  S1G3UT  is  calculated  for  each 
plate  comprising  the  panel,  assuming  the  panel  consists  only  of  that 
plate  thickness.  Thus  each  current  plate  thickness  within  the  trana« 
verse  panel  is  tested  far  buckling.  It  S1GBK3?  is  less  than  S1GBUT,  the 
plate  thicknesses  remain  as  the  current  thickness  1(1,1,).  If  S1GBKT 
is  greater  than  S1GBUT,  the  plate  thickness  is  increased  to  the  mean 
value  between  the  current  thickness  T(X,L),  and  the  required  thickness 
TL(J),  the  thickness  required  to  withstand  S1G3KT.  When  S1GBKT  is 
within  100  p.s.i.  of  S1GBUL,,  the  panel  satisfies  buckling. 

(BUK. )  is  executed  at  three  points  in  the  program:  at  the 
beginning  of  Cycle  B,  and  at  the  end  of  Cycles  A  and  B.  At  the 
beglnaaag  of  Cycle  B,  (BUK.)  determines  which  plate  panels  buckle*  so 
that  the  linear  program  in  (CYC.)  knows  which  panel  thicknesses  to 
change .  (BUK.)  is  executed  at  the  end  ©f  Cycles  A  and  B  to  determine 
if  any  plates  buckle  after  being  changed  by  the  linear  program. 

Four  possible  situations  exist  when  (BUK. )  ie  executed  at 
the  end  of  Cycle  A  and  B. 

X»  Limiting  stress  satisfied  and  no  plates  buckle. 

31*  Limiting  stress  not  satisfied  and  no  plates  buckle. 

2X1.  Limiting  stress  not  satisfied  and  plates  buckle. 

IV.  Limiting  stress  satisfied  and  plates  buckle. 
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3En  cases  I  and  XI,   tfce  thicknesses  entering  (BOX.)  a«T<? 
returned  to  (CYC.}.  In  case  III,  (BUK. )  increases  tbe  thicknesses  of 
the  plate®  that  buckle  in  the  manner  described  previously,  ajid  tbe 
thicknesses  are  returned  to  (CYC.)*  2»  case  XV,  (BUK.)  will  cycle 
until  all  plates  satisfy  tbe  buckling  stress.  Tbe  thlckae.  ses  of  tbe 
panels  tbat  buckle  are  increased  in  tbe  manner  previously  described. 
tbe  new  eress« sectional  area  properties  of  longitudinal  coz»tl&uou& 
members  &re  calculated  by  (AMIPLT.)  aad  (AM1SXF. ),  **&&  ^ae  »''    ""«•» 
distribution  is  determined  £>y  (9XRES1.).  Bew  stiff ener  a     tags  are 
calculated  by  (STFHKR.)  aad  then  {VHKHBS>)  is  executed  ta  ke&p  tbe 
plate  thicknesses  up«te«date  wltb  tbe  new  stress  sc     s.  Tbe  new 
stress  distribution  and  pansl  thicknesses  are  returned  to  tbe 
beginning  of  (BUK.)*  where  recycling  occurs.  Only  after  ©11  plate 
thlckeasses  meet  tbe  instability  requirseents  will  cycling  cease,  and 
control  ratunuB  to  (CYC.)*  Far  all  cases,  when  tbe  p&a^l  thleknesses 
are  returned  te  (CYC),  they  are  printed  out  under  tbe  beadl&gs: 

"BBggHQ  G0VEES3"  for  Cycle  A  and 

"BOCKLXKG  QOVERSS  for  Cycle  B» 
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LIST  OF 
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BT 
CHG<I) 

cow 
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SJLGBUL 
SlCSUT 

TB 
TL{I) 

TOP(K} 

TCTCHG 
TP 


late  panel  In 

sec  between  the  ac1 
acting  on  a  panel , 


Dimension  of  a  plate  panel  perpendicular  t©  the  l*>adv 

Tbe  number  of  a  plate  which  Is  located  below  the 
neutral  axis  and  buckles. 

The  total  number  of  plates  which  are  located  below  the 
neutral  axis  and  buckle. 

The  number  of  a  plate  that  buckles . 

The  constant  "Kn  used  in  &Jfyanss  Formula  which  is  a 
function  of  the  edge  restraint  conditions  and  the  boundary 
conditions . 

Applied  primary  stress  acting  on  a  longitudinally  framed 
panel o 

Critical  buckling  stress  determined  from  Bryan 3e  Formula* 

Critical  bu-kling  stress  determined  from  Montgomery  s 
Formula. 

Actual  thickness  ©f  a  longitudinally  framed  panel. 

Required  thickness  of  a  longitudinally  ©r  transversely 
framed  panel  to  withstand  '-'^viprcselve  loading. 

The  number  of  a  plate  wfcich  is  located  above  the  neutral 
axis  and  buckles. 

Total  number  of  plates  that  buckle. 

The  total  number  of  plates  which  are  located  above  the 
neutral  axis  and  buckle. 
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